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ABSTRACT 


This thesis describes a Bayesian method to determine 
the number of samples needed to estimate a proportion or 
probability with 95% confidence when prior bounds are 
placed on that proportion. It uses the Uniform [a,b] 
distribution as the prior, and develops a computer program 
and tables to find the sample size. Tables and examples 
are also given to compare these results with other 
approaches for finding sample size. The improvement that 


can be obtained with this method is fewer samples, and 


consequently less cost in Weapons Testing is required to 
meet a desired confidence size for a proportion or 
probability. 
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I. INTRODUCTION 


PProbabimemey, 1Smerelative, In part to 22. ignorance; in 
part to knowledge." (Ref. l:p. 140] 

This is the epitome of kLaplace's interpretation of 
probability, stated in the 1951 translation of his book A 
Philosophical Essay On Probabilities. Thieme ODplecsof this 
thesis is to estimate a probability. In particular, we 
will try to answer the question of how many trials are 
necessary or what should be the sample size to estimate a 
MeecPpertion Or probability from a set of Bernoulli trials. 

In many forms of Weapon System testing, sampling is not 
done sequentially, and the number of items to be tested 
must be specified before testing begins. Clearly, enough 
weapon systems or components must be tested to furnish 
reasonable confidence in the resulting estimate of, Say, 
system reliability. On the other hand, since testing is 
expensive and often destructive (e.g., missile launches), 
the sample size should be no larger than necessary. 

Many measures of effectiveness for military systems are 
mame the form of proportions, or probabilities of an 
attribute occurring. Some examples are 
System Reliability, 

Hate Propabi lity ; 
Launch Probability, 


Detection Probability, and 
. Fraction Defective. 


Ul im Ww DO F 


In such cases, testing may often be described as performing 
a set of independent Bernoulli trials. 

The problem is stated as follows: how many Bernoulli 
trials must we conduct, so that with a certain level of 
confidence, we can estimate the appropriate proportion or 
probability. A way to approach the problem is given by the 
definition of a confidence interval: 

A confidence interval for an unknown parameter gives an 
Indication .vor the numerical value of our unknown 
parameter as well as a measure of how confident we are of 
that numerical value [Refs 22p. 3255— 

Given a desired confidence interval size tcpaia 
proportion or probability, we wish to know the numberggs 
Samples needed to provide a confidence interval of that 
Size, and in this thesis we will produce tables and a 
computer program to assist in finding that sample size. 

We will discuss two methods for the above calculations 
and we will compare the results. The first and well known 
one from classical statistics bases the estimate upon a 
Simple random sample, and confidence intervals and sample 
size are explained in the next chapter. The second method, 
and primary focus of our study, is the Bayesian one. The 
basic advantage of this method is that it makes better use 
of the existing experience of the experimenter and his 
knowledge of the phenomenon being studied. It aggregates 
the information prior to the execution of the experiment 


with the observations after. This different concept uses 


Bayes' Theorem, and may result in smaller sample sizes 
while providing the same sized confidence interval. 

In Chapter III, we will describe Bayes' Theorem with 
Eheompntonm, sampling, and posterior distributions, will 
explain the use of the experimenter's prior bounds on the 
proportion and the choice of the Uniform distribution as 
prior, and will give the derivation of the posterior 
distribution and its properties. Then, in Chapter IV, we 
will calculate the sample size needed to estimate a 
proportion and we will compare the results with the 
classical method. We will explain the computer program 
used for the Bayesian results and will provide tables and 
examples to assist the reader. The final chapter will 
Summarize our work, and suggest additional applications of 


Bayes' Theorem to reduce the cost of weapon system testing. 


fee SAMPLE SIZE TO ESTIMATE A PROPORTION 
USING THE CLASSICAL METHOD 

In this chapter, we will explain the classical method 
to find the sample size to estimate a proportion. First we 
will find a point estimate of our proportion or probability 
which is an estimate given by a single number. Then we 
will find an iauinterval estimate, given by two numbers 
between which our proportion must be considered to lie. 
Interval estimates provide an indication of the precision 
or accuracy of an estimate and are therefore preferable to 
point estimates. Finally, we will use this confidence 
interval to determine the number of samples needed tc 


achieve a particular interval size. 


ye THE POINT ESTIMATE FOR A PROPORTION 

Generally, an estimation problem consists of the 
Manipulations we might make of the observed values ina 
Sample to get a good guess, or estimate of the value of an 
unknown parameter or parameters. 

In our case, we have a sample of n items. The 
probability of occurrence of an event (detect a defective 
item), called its success, is p while the probability of 
non-occurrence of the event is 1- p. We inspect all the 


nm items and count the number of successeS aS a Sequence 


Cmmmece pendent = bernoulli trials. Let xi be the outcome of 
each trial, where 


1 if we have a success, and 
0 if otherwise, 


xi 
xi 


and let x be the total number of successes. Then the point 


eetimate for our proportion will be the sample proportion 


Nod 
a 


415i __ = (2.1) 


MH 
tl 
tl 


where x follows a binomial distribution. The distribution 
of sample proportions has mean Wp and standard deviation op 


given by 


Lp = p and op = p(li-p) (242) 


feet. 3:p. 142). 

For large values of nn, the distribution of sample 
proportions iS approximately normally distributed. itigl 
marticular: 


The normal curve gives an excellent approximation to the 
binomial distribution when p is close to 0.5. EDwtace, 
for p = 0.5, the approximation is good for n as small as 
10. As p deviates from 0.5, the approximation gets worse 
and worse. On the other hand, for values of p 
picdnaticGantiy G@aitferent from 0.5, “the approximation*of 
the normal distribution to the binomial distribution gets 
better, the larger the value of n. Even if pis as low 
Foe O.10 or as high as 0.90, if n runs above 50, the 
normal approximation does not give bad results. Below 
0.10 or above 0.90, the Poisson distribution is commonly 
used to approximate the binomial distribution, although 
the normal distribution still does fairly well so long as 
enee > [Ref. 4:p. 100]. 


B. THE CONFIDENCE INTERVAL FOR A PROPORTION 


Let os be the standard deviation of the sampling 


distribution “of 9a = statist se S- ise the sampling 
distribution 1S approximately normal, we can expect to 
find, or we can be confident of finding, an actual sample 


Statistic S lying in the interval E[S] - 30s to E[S] + 30s 
about 99.73% of the times. Because of this we call this 
interval the 99.73% confidence interval for estimating 
E{S]. The end values (S + 30s) are the confidence limits. 
Similarly, S + 1.960; and S + 2.58 os are 95% and 99% 
Eontidence limits Lor Ws. The percentage confidence is 
called confidence level and the numbers 1.96, 2.58, etc..., 


in the confidence limits are called confidence coefficients 


and are denoted by 2c. For this study, we will work with 
the 95% confidence level, the normal approximation to the 
Binomial, and the corresponding 1.96 confidence 
coefficient. 

If the statistic S = Bp is the proportion of successes 


from a sample size n drawn from a binomial population in 
which the proportion or probability of success is p, the 


confidence limits for p are [Ref. 4:p. 572] 


Dat Zelaya) ‘ (2eF 
n 


We can compute the confidence limits of Equation 2.3, 
uSing the point estimate for our proportion from Equation 
2.1 and so the actual probability will lie in the interval 


6 


Pie le Oo Dil =e es PS Pee 1.96 pil = py (ee } 
Nn n 


with a 935% confidence level. For “example; 2£ from a 
population we inspect 30 items and 6 are found defective, 
we can say that we are 95% sure that the true value of the 
defective proportion p will lie in the interval calculated 
from the above Equation, where 6 = 6/30 = 0.2, 
Cre On 4 <0. 2045,0.14, 
OG O206. << p <= 0234. 
hie mm inececavalamsilze 12S Om34°"—- O0.06,= 0.28, and this 


becomes smaller when the sample size increases. 


C. DETERMINING THE SAMPLE SIZE FROM CONFIDENCE INTERVAL 

Let's state our problem again, as it was discussed in 
@hnapter I. How many items must we test so that with a 
certain level of confidence, we can report the reliability 
of this type of item. The certain level of confidence will 
be 95% for this study. 

One measure of the effectiveness of a sampling effort 
1s the accuracy of the resulting estimates. In our case of 
estimating a proportion, accuracy is reflected by the size 


of the resulting 95% confidence interval, or 





If the experimenter is willing to specify the size of the 


confidence interval on p that results from his testing, 


then his requirement may serve as a basis for specifying 
sample size. 
Let 2A be the desired 95% confidence interval size. 


Then our proportion will lie between 
Db — AS plsee ee 
and the interval size is 6 - A to 6 + A or 
Bp tA 
From Equation 2.4 we have 


(i Diy 22) , {2 oe 


and the sample size n can be determined by solving Equation 


20 (Rei. Wiap. 247 ) 


f 


_ (4,28) " =. (2.6) 
umes 


The sample size increases with the accuracy that we want 
for our estimate. Better accuracy means smaller interval 
Size, that is, smaller A and thus from Equation 2.6 bigger 
sample size. Also, the sample size is proportional to the 


square of the confidence coefficient, which reflects the 


desired confidence level. Finally, the sample size depends 
on our guess for the proportion p, before we actually 
sample from the population. We find the first derivative 


en EGuatwoen, 2.6 to be 


dnvdpe- (1.967A)2 (1-2p) 


and the second 


G-taei- eae 1. 967A) 2 p 


This is negative; so the value of p that makes the first 
derivative zero maximizes n. This happens for p = 0.5. 
Thus, our worst case where we need the maximum sample size, 
is found to be when we guess that half of the population is 
defective or that we have 50% chance to detect a defective 
item. In this case we need to sample n = (1.96/A)2 (0.5) 


mo. 5) or 


n = 0.9604/A?2 .on 7} 


items. The sample size decreases when the probability of 
Success increases from 0.5 to bigger values. Finally, this 
Gives an interpretation of the requirement: the value of 
2A is the largest confidence interval that the experimenter 
is willing to have result from his sampling. The value of 
h given by Equation 2.7 will guarantee that his requirement 
is met. 

Table 1 shows the required number of samples to obtain 
different 95% confidence interval Sizes for various 


reliabilities. 


TABLE 1. NUMBER OF SAMPLES TO OBTAIN 95% 
CONFIDENCE INTERVAL SIZE 
(Results rounded up) 


Interval Probability Of Success = p 

Size = 2A O25 O26 Oe 7 Ge Oe, 02995 
O205 leo i, 2 one 2 oe 984 554 150 
O16 535 367 S22 246 139 38 
G15 Va 164 144 110 62 ai 
0.20 oF 93 81 62 35 1G, 
0.25 62 60 52 40 Pies 6 
OP 0. 43 41 36 28 16 5 


| 


From the above table, we see that if we think that our 
probability of success will be p = 0.8 and we want to be 
+0.10 accurate with 95% confidence level, we have to sample 
62/1tems- 

The numbers of the above table are used to construct 
the <“qraph Gn Figure ae Here, we visualize better the 
previous discussion about the changes of the sample size 
because of interval size and probability of success. 

In Chapter III, we will solve our problem with the use 
of Bayesian methods, which give better results, that is, 


smaller numbers of samples than those in Table 1. 
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III. THE BAYESIAN METHOD AND ESTIMATORS 


In this chapter, we will explain a Bayesian method to 
find the sample size to estimate a proportion. To do this, 
we will first recall Bayes' Theorem and we will use it in 
our problem. Then we will explain analytically the three 
parts) of the Bayesian result we found: the prior 
distribution, the sampling, and the posterior one. We will 
state our reasons for the selection of the Uniform in the 
interval [a,b] as prior and Binomial as sampling 
distrabiut vour After that, we will derive the posterior 
distribution and its first two moments and we will find the 
Bayes estimators. Finally, we will explain the assumption 
we made in order to use the posterior distribution to 
calculate the 95% confidence interval of cur properterzem 

Inferentially, the Bayesian method permits the use of 
the knowledge and past experience of the experimenter, 
before observations are taken. Those, in combination with 
the sampling results, may give a smaller number of samples 
to estimate a proportion than that given by the classical 


method. 


A. BAYES' THEOREM 

One different method to estimate a proportion is to use 
Bayes' Theorem. Let us explain the procedure stating 
Bayes' Theorem first. 


ee 


SUimcosce = that —=ttyuene,-.-.., An are mutually exclusive 
events whose union is the sample space S, i.e., they form a 
partition of the S and one of them must occur. Then if A 


is any event of S, we have the following Bayes' Theorem: 


P(Axr) P(A/Ak) 


P(Axk/A) = a"; Sr ” (Sea 
= P(Axr) P(A/Ak ) 
k=1 
Consider now our problem. If alot has a defective 
proportion p, then the probability that a sample of size n 


will contain exactly X defective items is, for relatively 


imieoe 1OtS, approximately [Ref. 4:p. 558], 


ee a Pk = p)a>* Be) 
Suppose now that p is itself a continuous random variable 


with density function f(p), where 
\. Pi poadp = 1 


Mten the jOint probability that for a given lot, (1) p will 
fall in the interval p to p+ dp and, (2) that a sample 
size n taken from this lot contains X defective items, is 
Ene product 

P(X,p) = P(X/p) £(p) : 
According to Bayes' Theorem above, for the continuous case, 
the probability that the p that produced the given X lies 


ja the interval p to p + dp is 


PAA PD) &(p) 
Pp/ A, (os) 


|, Pix Dp) £ip) ap 


re 


The denswty Bimct tonne, is called the prior 
probability distribution and the probability P(p/X) is the 
posterior probability distribution. The third part of the 
above Equation 3.3, the sampling distribution P(X/p), is 
the probability function from which we will take the X 
items. Because we count successes in repeated n Bernoulli 


trials, this is Binomial as inv Bawecton. s7-- 


B. SELECTION OF THE PRIOR Dis tk tsuneo) 

The prior distribution of a parameter Pp 1S =a 
probability function oor probability expressing our degree 
of belief about the value of p, prior to observing a sample 
of a random variable xX whose distribution depends on p 
[REE 2 Sp) ooo. In other words, we can assign a prior 
GLsSte i bute von to a parameter Pp when we have’ enough 
information about the relative frequencies with which p has 
taken each of its possible values in the past. For 
example, suppose that the proportion p of defective items 
in a certain lot is unknown. Suppose also that this lot is 
made from a manufacturer who has produced many such lots in 
the past and that detailed records have been kept about the 
defective fractions in these Ors. The relative 
frequencies for these past lots can be used to estimate a 
prior distribution for p, which can be used in our Cereaus 
Pot: 

Different distribution functions can be characterized 
as “priors”. As examples, for a bounded variable p, we 
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mMentionmene Uniform distribution on the interval [0,1]. a 
triangular shaped distribution, and the Beta distribution 
with various parameter values. The Beta distribution for 0 
$p s 1 was used as the prior in Ref. 6, where the sample 
size problem for a Bayesian confidence interval was also 
addressed. 

On the other hand, we must note that the prior 
Srmseributton 1S a Subjective probability distribution in 
the sense that it represents an individual experimenter's 
information and subjective beliefs about where the true 
waive Of =p 15 likely to lie". [Ref. 5:p. 314] Often the 
best prior information about the parameter p may simply be 
bounds on p, wherein the experimenter can only say that p 
will not exceed some value b, and will not be less than 
some value a. The density function that is reasonable to 
combine with experience expressed as bounds on the unknown 
parameter seems to be the Uniform distribution on the 
interval [a,b] since it "distributes our ignorance equally" 
in the prior known interval [Ref. 4:p. 560]. 

The Uniform density function and prior distribution for 


this study is the Uniform [a,b] 


Rona “Sp Sb 
f1 (p) = (34) 
0 otherwise, 


where 
Os as ps b si 


ie 


Note that the Uniform [0,1] distribution belongs yea 
class of Beta (r,s) distributions when both Beta parameters 
are Q. 

It is valuable to remember here that the Beta density 

function (in the form that we use extensively later) is 

(y+) s +2) 
oot rae — arelas for 0 s x =a 
Bi ams? (Ss tee) 
Lx) = 
0 otherwise, 


where 
r > =] anad's) > =2 


C. DERIVATION OF THE POSTERIOR 

The posterior density as it is expressed in Equation 
3.3, Bayes' Theorem, iS simply the conditional density of 
p, given the sample values. It "expresses our degree of 
belief of the location of p, given the results of the 
Sample’. (Ret. 2p. “5561 

To derive our posterior distribution f2(p/x), we 
substitute into Equation 3.3 the prior as the Unitomm 
density from Equation 3.4, and the sampling distribution as 


the Binomial from Equation 3.2: 


2 a) 
oy (1 _ p)n-*x —_—_—_—_——_—_ 
x ) b a 


fo (p/x) cr ; 
n x —_ n-x — + = 
. ~ DucGe p) ; dp 
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‘ 
9) 


where 


IA 
be 


0 < a 


lA 


Pps b 
If we cancel out terms, we have 


p* C4. — p)a-x 
fe (p/x) = 


b 
I pt) -1p)"-< ap 


If we multiply numerator and denominator by the same 


number, we have 


r (n+2) 
ee p* (lp) ws ~ 
hoot er (nxt) 
TESS LO, DG 
x r(nt+2) 
| ae 2. ays XCD 


*  T(x+1) I (n-x+1) 


and we notice that the denominator is the area under the 


curve of a Beta distributed random variable with parameters 


ma x + oe 
and San - xX + 1 F 
over the interval ato b. Thus, we have 
Tr (n+2) 


a ee en TX fay) 8 = x 
r(x+1) I(n-x+1) 
eS ee ene ; 
Fa (b) - Fs (a) 
where Fs is the CDF of a Beta (r,s) distribution. 
We also note that the numerator has the form of a Beta 


density function with the same parameters (r,s) but the 


argument p is defined to be ecm es © 


i 


So finally, our posterior @istribution becomes 


Egutpy 
{2 (p/x) = ———— eee (3259 
F3 (0) een ae 
where O S$ a2 SDs ee 


F3 : CDF of Beta (r,s) 


and f3 (p) has the form of Beta (r,s) , 
where ro re 
and S =) — >. 
If we Lee use = 1/([Fs (b) - Fg3(a)], this posterien 


distribution has the functional form of a Beta density 
function with parameters (r,s), multiplied by a positive 
constant c 21.0, for a random variable p bounded by the 
boun@s (a,b) of the Uniform prior @dastrimeucio. 

Let us illustrate with an example uSing the above 
Conclusions. Suppose that an experimenter uses past data 
and puts bounds a = 0.2 and b = 0.8 on the probability that 
a defective item is located in a sample of n = 10 items. 
His prior distribution is Uniform for arandom VvVariableme 
bounded between 0.2 and 0.8. After the inspection of all 
the items, he counts x = 5 defective. Then from Equation 


3.5, we have r = 5 + 1 = 6 ands = 10 - 5 +1 = 6. We also 


have from Beta CDF that F3(0.8) = 0.98834, F3(0.2) = 
0.01165 and c = I1/[F3(0.8) - F3(0.2)] = 1/05397G63— 
1.02386. Then, his posterior distribution, 1S a )foernmeen 
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Berane Oo) mubtipided by 1.02386, for a random variable p 


bounded again between 0.2 and 0.8. 


D. THE POSTERIOR DISTRIBUTION AND BAYES' ESTIMATORS 

We have shown above an example of the density function 
Of the posterior distribution, for specified Uniform prior 
and Binomial sampling distributions. We will now calculate 
the mean and the variance of the posterior distribution in 
the general case of Equation 3.5. Let c be the constant 
factor 1/[Fs (b)-Fs (a)] in the Equation 3.5. Then 

i tee) 
E{p|x] = \. be nnn (Lp) ~ x dp 
Mex hl eae 1) 

If we combine terms and multiply numerator and denominator 


by n+ 2 and x + 1 respectively, we have 








Sot < eet 3) 
E{p |x] = c | a a) x 1 ( L=pH) i> xX dp 
n+ 2 “T(x + 2) [(n - x + 1) 
> ae naa b 
a | fq (p) dp t 
n+ 2 : 
Where £4 a form of Beta (r = x +2, s=n- xX + 1) 


and asp sb. 


PMbostituting for €, we have 


x +1 Fa-(b) — FPas(a) 





E[p|x] = 7 (2mey) 
n+ 2 F3 (b) - F3 (a) 

Where Fa : CDF of Beta (r = x+2, s =n- xX +1) 

and F3 : CDF of Beta (r =x +1, s=n-xtti). 


des, 


It is worthwhile to give a different form of Ehe mean 


From Equation 3.6, we have 





bo TT (ae a 
x+1je Pix + 2) fm =e) ) ee dp 
E{p|x] = 
n+ 2 [Pe r(n + 2) ; ~~ 
I, r(x + 1) Fin - x + 1) Pp (1—p) x dp 


Since the arguments of the Gamma functions are integers, we 
can substitute with factorials and pull them out of the 


integrals, giving 





pee ea ? are “pete 
is os (xc 1) tener |. p**+ (1-p) dp 
Be ple 
+ 2 - b 
° aera Nase a |, p* (l—pit sale 


Simplifing the constant terms, we have 


b 
| p* +! (1l-p)"°-*dp 
a 


(mn + 2)! x! (x + 1) (n -— x)! 

Bie alte 
b 
Pare OMe) (ee Ga sz) | | p* (1=p)"oap 

OF 

(nm +2)! te 1) eee |, py +) (Lopes * de 
E(p |x} = ————_— 

b 
in + 2)!) (x + 1) thee I. px (1-p)"-*dp 


which obviously gives 


b 

I. px¥t! (]-p)8-*dp 
b 

| p* (1l=p)">*dp 
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me Calculate the Variance, we use the well known result 
Var[{p] = E[p?] - E[p]*. Working as we did for the mean, we 
have that 


Xx +1 x + 2 








E{p? |x] = c (Fs (b) - Fs (a)) ; 


10 @ i eee 
where Fs is a Beta CDF with parameters r= x +3 and 
Seen 6 -CUCUCU KC CU. Then the variance of the posterior 
distribution is 


Sct ae 

















Var(p|x) = c (Fs (b) - Fs (a)) 
Queer, 2 ay oer Ss! 
pie ap IL 2 
= ¢ (Fa (b) - Fa (a) ) ; 
Fe ee 
or 
x +1 x+ 2 ->Fs (b) —- Fs (a) 
Var(p|x) = SSS 


n+2 n+3 #Fs(b) —- Fs (a) 


Meee ae) o— haat a) \ < 





(Brie) 
nectar Fs (by > F2 (a) 


The commonly used point estimate 6 for a proportion in 


the Bayesian method is the mean of the posterior 
Gaestribution. However, if the posterior density for pis 
not symmetric, other measures of the middle of the 
posterior might also be used as the point estimate. Two 
such measures are the mode of the posterior (which 


maximizes the posterior density) or the median (which is 


eu 


the value that splits equally the area under the density 
Gurvel = 

As in the classical method, interval estimates are 
preferable to point estimates. In this Bayesian method, it 
1s easy to construct interval estimates of the proportion 
Da A 95% confidence interval is provided by the 2.5 and 
97.5 quantiles of our posterior distribution of Equation 
3.5. Thus, the interval estimates depend on the subjective 
bounds of the prior Uniform distribution (a and b), the 
Sample size ns, and the number x of successes from the 
sampling (Binomial distribution for sequential Bernoulli 
Pera ls): Before sampling, we know a and b, but we do not 
know n and x. Note that our problem is how big should be 
the sample size n. 

This results in a need to guess the number x of 
Successes before we actually sample. We recall from the 
definition of the mean of a random variable, that this 
number would locate the center of gravity of the 
distribution of the random variable and thus, is a likely 
candidate if we have to give a single number as our guess 
of the value of the random variable. So, prior to 
sampling, a “good" guess for x, the number of successes, 15S 
the mean of the prior Uniform distribution multiplied by 


the sample size n, or 
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MecinGmeetiits, OULD MDeSterior distribution in Equation 3.5 
becomes, for purposes of determining sample size n, 
f£3 (p) 
Poa 22) = i (so. 37 
F3 (b) - F3 (a) 


where O's avsep s Db Ss 1 ; 


F3; : CDF of Beta (r*,s*) 


and f3;(p) has the form of Beta (r*,s*) F 
at b 
where Ty FD 
2 
aot. D 
and ze aes 9 a eae Oo a oe A 
2 
Momeeduation 3.8, the posterior distribution has the 


functional form of a Beta density function with parameters 
fer, S*), multiplied by a positive constant ce2 1.0, for a 
random variable p bounded by the bounds (a,b) of the prior 
Wmrtorm distribution. 

BEOnmmnowmrcon Hmcduation 3.8 will be our posterior 
distribution which we use in the next chapter as the base 
to develop a procedure to calculate the interval estimates 
and sample size. In the next chapter also, we will discuss 
the computer programs that we used to find the sample sizes 


to estimate the proportion p. 
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LV: SAMPLE SIZE TO ESTIMATE A PROPOR 
USING THE BAYESTAN METHOD 

In this chapter, we will explain how to Use se 
Bayesian method with a Uniform prior in order to find the 
sample size to estimate a proportion. The experimenter has 
some prior information about the unknown proportion in the 
form of upper and lower bounds on the unknown proportion. 
We use them as bounds for a prior Uniform distribution and 
we wish to determine the sample size he needs, based on the 
accuracy he likes. 

First, we will derive the bounds of the Bayesian 
confidence interval, where the proportion to be estimated 
should lie with 95% confidence level. This interval leads 
to the necessary sample size. Then, we will discuss the 
computer programs we have used in this procedure. Finally, 
we will provide tables and examples to assist the user to 
find the sample size that meets his goals, and to visualize 
the advantage of the Bayesian method in giving smaller 


samples than the classical one. 


A. THE BAYESIAN CONFIDENCE INTERVAL 

Once we have obtained our posterior distribution, we 
can construct an interval which contains 100(1-a)%* of the 
posterior probability. Our posterior distribution, given 


by Equation 3.8, has the Beta (r*,s*) form, multiplied bye 
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constant for arandom variable p bounded between a and b. 
We let a be 0.05 for this study, and thus want a 95% 


confidence interval size. 


Let p.lo and p.up be the lower and upper bounds of the 
confidence interval. Letting the area between a and p.lo, 
the lower bounds, be a/2 = 0.025 of the whole area under 


our posterior density function, fz2(p/x), we have 
Fo (p.lo) = a/2 = 0.025 


From Equation 3.8, we have also that 


p lo 
Fo (p.lo) = [i fe (p/x) dp, 
or 
i p.lo 
me ——— | 22 Gon cl, 
F3 (b) —- Fs (a) 7 


where F3(p) is Beta (r*,s*). 


Thus, we have the equation 


1 p.tlo 
| f3 (p)dp = 0.025, 
Fs (b) - Fs (a) 
Long 
Fs (p.lo) ~ Fa (a) = 0.025 [Fs (b) - Fs (a)], 
which finally gives 
Pompeo) = 0.025 Fs (b) + 0.975 Fs (a) x Ca) 
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Similarly, letting a/2 = 0.025 be the dareawindenwaaa. 
posterior density function fz2(p/x) between p.up and b, we 


have 


Fs (p.up) = 0.0975 F3 (Bb) +> 0.6258). (cee (4.2) 


From Equations 4.1 and 4.2, we calculate p.lo and oem 
and then by subtracting p.lo from p.up, we get the 935% 
confidence interval size. This is the Bayesian interval 
where our proportion to be estimated should lie 95% of the 
time. As we did with classical method, let us call the 
Size of this interval 2A; it is a measure of estimation 
accuracy... We will see that the sample size depends upon 
the prior bounds (a,b) and upon the interval size 22h. The 
decision maker uses his past information to state the 
bounds and his preference in accuracy to state the interval 
Size. The Bayesian method of this study gives the interval 
p.lo to p.up, where the proportion p lies with 95% 
confidence level. 

In the next section, we will explain how these values 


may be used to find the required sample size. 


B. DETERMINING THE SAMPLE SIZE FROM THE BAYESIAN INTERVAL 
In Chapter III, we derived our posterior density 
function which has the form of a Beta diStributicn wee 


parameters (r*,s*), multiplied by a constaat ecw 


ONS. 


f3 (p) 
1a SC) SS 
F3 (b) - Fs (a) 
where GOA< a < pesebys. lL. 
Fos eDF “or Beta “(* ,s*) 
and fs (p) has the form of Beta 
(lL SW: 
In Chapter III, we also explained why, for purposes of 


determining the sample size n, the parameters r* and s* 


take the values 


a+b 
r* = ig) aces ; 
2 
and 
a+b 
Sos i) ae de ; 
2s 


where a and b are the bounds of the prior Uniform 
feestribution. 

Once we have obtained the parameters of our posterior 
Sersceripution, UWSsing Equations 4.1 and 4.2 we need to 
compute the inverse cumulative distribution function at 
woes and 0.975 for a Beta with parameters r* and s*. This 
Will result in the lower and upper bounds of the 95% 
confidence interval. Then, if we subtract the lower from 
the upper bound, we can determine the size of the desired 
confidence interval. 

The above procedure is used in the APL program SAMPLE 


Hecated in Appendix A. This program computes the sample 


Pad | 


Size needed to obtain a 95% confidence interval where the 
probability or  preportionssscnould meee The program is 
interactive and requires the user to input the bounds of 
the prior Uniform distribution and the desired confidence 
interval size. Then it calculates the parameters of the 
posterior distribution and computes the confidence 
interval that is provided when we sample 10 items. 
Consequently, uSing a loop, it increases the sample size 
until the desired confidence interval size is reached. 
Finally, it prints the 95% confidence interval bounds and 
the sample size needed to obtain the required confidence 
interval. 

The program SAMPLE uses the subroutines BQUAN, NQUAN, 
and BETA located in Appendix B. These are APL programs 
designed at Naval Postgraduate School to compute the 
inverse cumulative Q2istribuGren function of Beta 
da Scrip tion. It must be noted that BQUAN often cannot 


compute the inverse cumulative distribution function for 


large Beta parameters. In our case, large parameter values 
mean large sample size. Thus SAMPLE was written to 
terminate its calculations at sample size 150: results can 


always be obtained for sample sizes at or below this value. 
This number can be increased, but in general SAMPLE cannot 


evaluate sample sizes greater than say 200. 
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C. TABLES FOR FINDING THE SAMPLE SIZE 

In this section, we provide a table to assist the user 
in finding the sample size that reflects his past knowledge 
with the prior bounds, and his preference for accuracy with 
the 95% confidence interval size. 

Table 2 was constructed by executing the APL program 
SAMPLE repeatedly for selected values of a and b, covering 
the whole range from 0 to1.0. Also, convenient interval 
sizes were used. 

The use of Table 2 is Simple. For example, suppose the 
user puts prior bounds 0.5 to 0.8 and wants the interval 
Size 2A to be 0.25. He looks at the part of the table for 
b = 0.8 and he finds the entry in row a = 0.5 and column CI 
Size 2A = 0.25. He has to sample 40 items to be 95% 
confident that the proportion p will be in a confidence 
iueerval of size 0.25. 

Before giving more examples of the use of this table, 


moms well to note that some entries in Table 2 are blank. 


One reason for missing entries 18, aS mentioned in the 
previous section, that SAMPLE generally can not evaluate 
Sample sizes greater than 200. Another problem in 


constructing Table 2 occurred with the APL program NQUAN. 
Its execution stops when the sample size is big (greater 
than 100) and the sum of the prior bounds (a + b) is 
fietween 0.7 and 1.3. Finally, another reason for blank 


meres Psethat 2A must be less than b - a. 
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TABLE 2. NUMBER OF SAMPLES TO OBTAIN 95% CONFIDENCE 
BAYESIAN INTERVAL 














5b =..05 b= be=aee 
Cl size ZA CI size 2A CI size 2A 
a O2025 0.05 O.20 715 0075 0.10 0.81 
0 617 299 82 249 140 40 
02025 336 275 oa 32 
a6, es 299 156 
02075 318 139 
| 
b= 0.3 b= 0.4 
CI size 2A CI size 2A 
Ot e702 15 O02 0.25 20m On) 0.15. 40nR 0.25 Ge 
196 87 42 2 24555 206 60 36 18 
205 Peo al 94 a] 118 64 35 is 
nak 244 93 LAAs: 66 30 
as, 26 
bo= 025 b = 0.6 
126 70 44 29 141. 78 49 29 
s05 134 a4 46 30 147 82 op 35 
Ek ai 78 48 29 is. 85 53 36 
re 153 82 40 89 55 35 
oe 90 44 
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(TABLE 2 CONTINUED) 








ae Oe), b = 0.8 

CI size 2A CI size 2A 
a Opt somes 0.3 OFf.25 0.3 Ce ieeeoeeros O20. 25 0.3 
0 53 85 53 36 89 56 38 
o.05 87 55 oe | 91 57 39 
7.1 89 56 2g 92 Sis eis) 
Ge a2 58 SPS, 168 93 59 40 
O23 168 93 58 36 92 58 39 
0.4 90 44 89 55 a5 
o.5 153 82 490 
o. 6 126 
a b= 0.9 b = 0.95 
0 Be 58 a, 93 58 40 
oe 05 93 58 40 93 59 40 
o-1 168 93 59 40 93 58 40 
Oo. 2 92 ok: 39 91 5 39 
G3 89 56 38 87 S\2) a7 
0.4 53 85 56 36 147 82 op 25 
G5 141 78 48 29 134 74 46 30 
0.6 126 66 30 118 64 25 15 
O. 7 244 93 PAPA 25 Si) 
Gr 8 1516 





on 





O82 


O72. 








O05 


OTS 





OF 


220 


209 


io 


ivegh 
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(TABLE 2 CONTINUED) 


b= 70.975 
CI size 2A 
0205 
495 
226 
b= 0,975 
CI size 2A 
0 LS 20 2 Sao 6) gal 
168 2 58 40 
oa 53 40 
93 58 40 
90 a7 Bh) 
52 86 54 7 
144 80 50 34 
130 1 45 PPS 288 
a3 62 36 17 245 
3) 41 8 196 
v2 cS oe 
a2 140 
ae 
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b= 1.0 


CI size 2A 


O05 


432 


200 


eZ 


b= ed 


CI size 2A 


0.15 0.2 0.2595mm 


168 93 
93 
a2 
89 

ee. 35 

141 78 

25 79 

108 60 

87 42 
gpl ay 
40 


og 


3 


58 


a6 


oe 


4g 


43 


2S 


12 


40 


40 


32 


38 


36 


29 


a2 
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We should mention also that all the programs used in 
Bouls study are ~writtem, 1n APL and can. be .run...on.any 
computer with APL capabilities. However, because of the 
extensive loops they have, they may require a significant 
amount of time. 

Let us illustrate with some additional examples the way 
that Table 2 can be used, andat the same time, let us 
compare the results it gives with those given in Table i 
tBeom classical statistics. 

The classical procedure requires the experimenter to 
state values for confidence interval size 2A and estimated 
probability of success p. The Bayesian procedure requires 
2A and prior bounds a, b. In order to be able to compare 
the results obtained from these two methods, we recall an 
argument in which our Bayesian procedure was based. The 
experimenter, using his past knowledge, states the bounds 
@eeene prior Uniform distribution. 

For any finite sample size, the Bayesian estimate is 

shaded" toward the prior mean, the best guess for © 

before any sample values were taken [Ref. 2:p. 566]. 
In our study, we use p in place of ®@. The mean of the 
prior is (a+  b)/2, the sum of the bounds divided by two. 
To compare the two methods, this number, from the Bayesian 
method, is used as the probability of success to enter 
Table 1 and find the suggested sample size from the 
classical method. Thus, if the bounds of the prior Uniform 


eee O0.5 to 0.9, we use the number (O253+ 0.9)72 = 0.7 as 
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the probability of success to enter in Table 1) andvenen 


compare the results. 


We proceed now with some examples to explain the use of 


the Bayesian Tables and how to find the sample size when we 


have the bounds of the prior Uniform distribution and the 


desired 95% confidence interval size. 


alee 


Example 1: Fraction Defective 


Suppose a lot of 10,000 items is received from a 
supplier; the lot contains p (unknown) defective 
items. Also, suppose that we have kept records for 
the past lots from this supplier andwe decide our 
Subjective bounds to be a = 0 and b = 0.4. How many 
items do we have to sample in order to be 95% 
confident for the fraction defective, with = Omen 
estimation accuracy? We look at the part of Table 2 
with b = 0.4, and we find the sample size of 60 items 
in the entry for row a = 0 and column CI size 2A= 
O24 If we look at Table 1, for probability of 
success 0.2 (column 0.8) and interval size 0.2, we 
find 62 items. The BayeSlan approach reduced a 
sample size by 2/62 or 3%. Note also that had we 
used the common textbook formula for sample size for 
a proportion, (Equation 2.7), the result is n = 97. 


Example 2: Hit Probabiliwve: 


Suppose the size of the load of a recently modified 
weapon system has to be decided; this system has 
(unknown) hit probability, p, against one of the 
targets it is designed for. Suppose also that we 
have data from the past firings with the old version 
and we decide our subjective bounds will be a = 0.7 
and b = 0.9. How many items do we have to fire in 
order to be 95% confident of the hit probability, 
with + 0.05 estimation accuracy? We look at the part 
of Table 2 for b = 0.9 and we find the number of 244 
items in the entry for row a = 0.7 and column CI size 
2A = 0.1. For the classical result, we look at Table 
1, for fp = 0.8 and 2A = 0.1, we find 246 items. With 
the Bayesian approach, we need 2 items less, i1.e., 
0.008% better results. Had we used Equation 2.7, the 
result would be 385. 


34 


SW Example S62 Detection Probabilrty 


Suppose the number of acoustic devices has to be 
defined in order to design a new type of sonobuoy. 
Each device has p (unknown) probability of detection. 
Suppose also that experience from similar hydrophones 
gives a probability between 0.1 and 0.5, and we like 
+ 0.15 accuracy. How many acoustic devices will we 
have to use in order to fulfill the requirements 95% 
of the time? We look at the part of Table 2 for b= 
0.5 and we find 29 devices in the entry for row a = 


0.1 and CI size 2A = 0.30. If we use Table 1 from 
classical statistics for probability of success 0.3 
(column 0.7), we find 36 devices. The Bayesian 
approach gives 7 devices less, that is 19.5% better 
results. Table 1 in column 0.5 (Equation 2.7) gives 
the number 43, almost twice as big as the Bayesian 
one. 


We have demonstrated above the use of the BayeSian 
tables ain finding sample sizes. The tables provide the 
most common bounds and interval sizes. If the user needs a 
Sample size for bounds and/or interval size that are not 
included in the tables, again the program SAMPLE can be 
used. The user interactively inputs his prior bounds and 


desired 95% confidence interval size and the output is the 


confidence interval bounds and the number of samples. An 
APL session, solving a problem in this case, is shown in 
Figure 2. Prior bounds are 0.40 to 0.75 and interval size 


is 0.25. The Bayesian method requires 54 items . 

In order to compare the answer with that from classical 
method, we use the Equation 2.6 for 6 = (0.40 + 0.75)/2 = 
0.575 and A = 0.25/2 = 0.125 and we find 61 items. [In this 
case, the Bayesian approach gave 7 items less or 7/61 = 


11.5% better results. 
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To assist the comparisons between the two methods, let 


us present the results of the examples in a table. 


SAMPLE 
ENTER A, LOWER BOUND OF UNIFORM PRIOR 
ble 

elt 


ENTER B, UPPER BOUND 
[eles 
e755 
ENTER 95 PERCENT C.L. INTERVAL SIZE (MUST BE LESS THAN B-A) 
ls 
25) 


95 PERCENT C.L. UPPER BOUNDT: 0.655320 24%20 
>> LOWER > QO.4467326643 
>> INTERVAL SIZE: Q.2491905779 


REQUIRED SAMPLE SIZE 3: 54 


Figure 2. An APL Session Using The Program SAMPLE 


TABLE 3. NUMBER OF SAMPLES TO OBTAIN 95% CONFIDENCE 
INTERVAL SIZE 


Number of Samples Improvement 
Example _a_ —b CI Size Bayesian Classical +e 
Example 3 0.1 0.5 0.3 29 36 19.5 
APL session 0.4 0.75 “0225 54 61 1145 
Example 1 0.0 0.4 0.2 60 62 3.0 
Example 2 0.7 0.9 Ont 244 246 0.008 
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We see that as the sample size gets larger, the *% 
improvement of the Bayesian method decreases. "The 
difference between the Bayesian values, and the classical 
approach, disappears as n increases" (Ref. 2:p. 573]. This 
happens because as the sample size gets larger, the 
posterior distribution becomes’ less dependent on the 
assumed prior and more on the sampling one. When values of 


mn are smaller, the Bayesian values may differ considerably 


from the classical. This situation underlines’ the 
importance of the prior distribution in that, for small 
sample sizes, the prior distribution must be chosen 
earefully. 


To complete the study of the Bayesian method, let us 


look again closely at Table Ze We see that the sample 
size n increases when the interval size 2A decreases. We 
see also, that for the same interval size and holding one 


bound fixed, the sample size increases as the other bound 
approaches 0.5. If we use the interpretation which we did 
before, i.e., to consider that the sum of bounds divided by 
two in Bayesian method is equivalent with the probability 
of success in the classical method, we conclude that 
Bayesian intervals behave exactly in the same way with 
enose in classical statistics. The discussion in Chapter 
II and its presentation with Figure 4, about the changes 
and the dependence of sample size in classical method, is 


also valid for the Bayesian one. 


ou, 


In the next chapter, we summarize our work and propose 
additional studies for the use of Bayesian methods to 
reduce the sample size and thus, the cost of weapon 


tLestinge 
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V. SUMMARY AND SUGGESTIONS FOR FURTHER STUDY 


In this chapter, we will summarize the procedure we 
used, working with the Bayesian method, to obtain sample 
sizes which are smaller than those given by classical 
Statistics. Finally, we make some recommendations for 
further research in using Bayesian methods to reduce the 
sample size needed to estimate a proportion or probability 


in any test field, and thus to reduce the cost. 


A. SUMMARY 

In this paper, we used a Bayesian method to obtain the 
number of samples needed to estimate a proportion or 
BPeonability. 

First, we described the classical method and explained 
the point and interval estimates. Using desired confidence 
interval sizes, we produced a table with sample sizes given 
from classical statistics for 95% confidence intervals. 

Then, we described the Bayes’ Theorem with the prior, 
Sampling and posterior distributions. We choose the 
Uniform [a,b] as prior. This, combined with the sampling 
Binomial, give a form of Beta distribution as posterior for 
a random variable bounded by [a,b]. We derived the 
Bayesian 95% confidence interval and produced a computer 
program to calculate the sample size. We provided a table 
and gave some examples to assist the user to determine how 
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to use these results to obtain smaller Sample .siz-ce 
Finally, we compared the results given by both methods for 
the same 95% confidence interval size. For small sample 
sizes, generally smaller than 100, the Bayesian method with 
the Uniform [a,b] prior improves the results and thus 
decreases the cost of tests based on sequential Bernoulli 
trials. 

Thus, when the decision maker has prior knowledge, and 
he wants to benefit from this, the Bayesian method of this 
study is recommended in order to reduce the number of 
items, and consequently, the cost. 

In the next section we suggest some additional studies 
based on Bayes' Theorem and on this paper, for even smaller 


sample sizes. 


B. SUGGESTIONS FOR FURTHER STUDY 
This paper uses the Uniform [a,b] distribution as 
Dr vor. This prior is easy to use, but it does not always 
give better answers than other prior distributions. The 
study by Manion approached the sample size question with a 
Beta prior distribution for a proportion bounded by —G=and 
a), For a quick comparison, we use an example. 
As an example, if the decision maker wanted the size of 
the 95% confidence interval to be 0.20 and his subjective 
bounds on the proportion were 0.14 to 0.86, the 
parameters on the Beta prior would be 4,4 and the number 
of observations needed would be 87 [Ref. 6:p. 42]. 
Our study with Uniform [0.14,0.86]” pricr 9giveou 


items. Possibly, it could be better if aeprior that 
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combines both studies and concepts were to be chosen, 1l1.e., 
a Beta for a prior bounded random variable. Another prior 
density function could be the triangular one, again for 
bounded random variable. The sensitivity of resulting 
sample size to the choice of bounds a and b could also be 
explored. 

An additional research task could be an effort to fill 
the blanks in the table of the Bayesian interval sample 
sizes of this study. This presupposes the development of a 
computer program that can compute the inverse cumulative 
density function of the Beta distribution for large 
parameters. 

Patd@weves dam addition tO this paper could be the 
development of tables for confidence intervals other than 
95%, such as 90%, 97.5%, and 99%. 

We hope that the chance to reduce the cost of sampling 
with smaller sample sizes to estimate a proportion, as 
meyer from this paper, will be beneficial to any authority 


dealing with tests of acceptance, reliability, etc. 
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APPENDIX A. THE APL PROGRAM "SAMPLE" USED TO COMPUTE 


C1) 

C2] 

[33 

C4) 

to) 

C6] 

(7) 

C8] 

C9] 

(10] 
Gala) 
[12] 
C13] 
C14] 
C15] 
C16] 
alga, 
C18] 
C19] 
[20] 
F244 
(22) 
[23] 
C24] 
C25] 
[26] 
b279 
[28] 
(29] 
(30] 
C31] 
F323 
(33] 
C34) 
(35] 
(36) 
(37] 
(38] 
[39] 
C4uo] 
C41] 
(429 
C43] 
Cuu] 
C45] 


SAMPLE SIZES FOR BAYESIAN INTERVALS WITH 
95% CONFIDENCE LEVEL. 


V SAMPLE [(Q] 

V SAMPLEsLEFsAT3N3C:D: LO; UP3E 
A THIS PROGRAM COMPUTES THE SAMPLE SIZE NEEDED TO OBTAIN A 
A BAYESIAN INTERVAL WITH 35 PERCENT CONFIDENCE LEVEL, BASED ON 
A A PRIOR UNIFORM (A,B] DISTRIBUTION. IT ASKS THE USER TO 
A INPUT THE PRIOR BOUNDS A AND B AND THE DESIRED INTERVAL SIZE. 
aA IT NEEDS THE APL PROGRAMS BQUAN, NQUAN AND BETA TO BE STORED. 
A IT TERMINATES ITS EXECUTION WHEN THE SAMPLE SIZE IS 2 150. FOR 
A BIGGER NUMBERS, THE VALUE OF N IN LINE 22 MUST BE INCREASED. 
A IF CONFIDENCE LEVEL DIFFERENT TRAN 95 PERCENT IS REQUIRED, LINES 
A 28 AND 29 MUST BE CHANGED ACCORDINGLY. 


"ENTER A, LOWER BOUND OF UNIFORM PRIOR' 

LF+O 

’ t 

"ENTER B, UPPER BOUND' 

RT+0 

t 

"ENTER 95 PERCENT C.L. INTERVAL SIZE (MUST BE LESS TRAN B-A)!' 
INT<0 


N<9Q 
CONT:?>FINx1N=151 
N«N+1 
C+14+Nx0.5x(CLF+RT ) 
D<N+2-C 

E«<c,D 


£O+(0.025xCE BETA RT))+025975x(CE BETASEED 
UP<«(0.975x(E BETA RT) )+0.025™ (CE SE A=) 
LO+E BQUAN LO 
UP+E BQUAN UP 


LOOP: *ENDx1((UP-LO )SINT ) 


+CONT 
END:' 
' ’ 
'95 PERCENT C.L. UPPER BOUND : ',®UP 
>> LOWER : ',3L0 
t t 
>> INTERVAL SIZE: ',®(UP-LO) 
t f 
‘REQUIRED SAMPLE SIZE : ',8N 
+0 


FIN: "SAMPLE SIZE IS GREATER THAN 150 AND EXECUTION TERMINATED! 
: 


t 


V 
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Et) 

(2) 

(3) 

C4) 

C5) 

C6) 

E77) 

(8) 

(9) 

(10) 
(11) 
(12) 
(13) 
C14) 
(15) 
(16) 
C17) 
(18) 
(19) 
[20] 
(21) 
22) 
(23) 
C24] 
[25] 
(26) 
[27] 
(28) 
(29) 
[30] 
31) 
[32] 
(33) 
C34) 
(35) 
C36) 
[37] 
C38] 
(39) 
C40) 
C41) 
C42) 
C43) 


APPENDIX B. THE APL PROGRAMS USED TO COMPUTE 
THE INVERSE CDF OF A BETA DISTRIBUTED 
RANDOM VARIABLE. 


Y BQUAN (0) 

UV VeA BQUAN PsE3:U3:S:D:b3:Z:DENS3:I3PP3sM3sX3F3:C23C33C4 

A IMPLEMENTATION OF CARTER, 1947, BIOMETRIKA FOR APPROXIMATE INVERSE BETA 
A 11/5/86 BEST FOR AC1)€2*AC2), AND SEEMS TO WORK FINE 

A 12/27/86 ADDED 2 NEWZQN-RAPHSON IZERATIONS: ADD MORE FOR GREATER ACC. 
+((L/A)<1)/SMALL 

E«NQUAN 1-P 

U<o"14+2xA 

S<4+/4U 

De-/4U 

Le ("3+Ex2 46 

Z+( (S42 )xEx (L424S )*0.5 )-Dx(L4+(546 )-S43 )-(Da2 )x ( (2485 )40.5)xEx (1148 *2 4144 
V<e414(4/0A )x*x2xZxJT<l 
LOOP: DENS+«A(1)]x(AC1J)!°> 14+4+/A)x(V*xAC1)-1)x(1-V)*A02)-1 

VeV-((A BETA V)-P)4DENS 

+((I«I+1)s2)/L00P 

+0 


Q 
A 
A MODIFIED 1/3/87 TQ USE MEAN AND STANDARD DEVIATION, AND NORMAL QUANTILE 
a 
A 


=o —_—<—< = Ts oe ee ee ee ae ee —_——_ —_———_ _——_ = = —_ee ewe = == 


SMALL: V«X<«(p,P)p0 

PP«+PSM«A(2)]4++/A 

XCPP/1pX1«((PP/P)*+(((AC1)<1),401321)/1,A01) )xAC1)$7~14+4+/A)*4A(1) 
XO (~PP)/1pXJ«1-((1-(~PP)/P)*+(((A02)<1),4A02]21)/1,A02) )xAC2]!$° 14+4+/A)*4#Al2) 
XO (X=1)/19X)]«1-1F 15 

+((f/A)21)/0NE 
START: F«(A(C1J]!~>14+4+/A)xAC1)= (X*A(01)-1)x(1-X)*A(2)-1 
C2<((1-AC1) )4+X )4+(A023-1)4+1-X 

C3«(2xC2*2 )4+((AC1)-1)4X*2)4(AL2)-1)4(1-X )*2 

Cu<(6xC2*3 )4+(7xC2x (C3-2xC2%2) )4+°((1-AlC1) )4X"3 4+(A02)-1)4(1-X ) 43 
Fe(P-(A BETA X))+F 

VeX4+F4+((C2xFx2 #2 )4+(C3x (F*3 46 )4+Cux (Fay )42u 

VO(V>1)/r1pV) +1 

+0 
ONE :M<1-M 

S« (Mx (1-M)4#144/A)*0.5 

+((A1F /A)=2)/44+0LC 

XCPP/19X]«M+SxNQUAN PP/P 

XO(XS0)/1pXJ]+«1E 15 

+START 

XO (~PP)/19X]4M4+SxNQUAN(*PP)/P 

XC(X21)/19X)]<«1-1F 15 

+START ; 

V 


43 


Cr) 
2) 
C3) 
Cu) 
ES] 
C6) 
7) 
C8] 
(9) 
(10) 
C11) 
fi2) 
(13) 
Ciu)j 
£15) 
(16) 
(17) 
C18] 
£19) 
[20] 
£21) 
£22] 


bid 
C28 
(3) 
Gee) 
tod 
(6) 
fad 
(8) 
ES) 
C10) 
Cit) 
[123 
(13) 
on 
£15] 
C16] 
7) 
C18] 
C19] 
(20) 


is 


V NQUAN [{0) 
V Z+NQUAN PsA: BsCsD:3sQ3sT73sS3R3F 


==] — eee ee ee ee — Se ee ee eet ee — oe eee Se ee — — = —— ss — en ee ee eet ee ee 


ESPECIALLY FQR EXTREME P VALUES, ADD QNE QR MORE NEWTON-RAPHSON LOOPS. 
+(v/(PsS0O),(P21))/ERR 

+(v/(C1Q«,P-0.5)s0.42))/34+05C 

S+Z+,Q 

+EXT 

T+(0.422|1Q)/Z+,Q 

+(F«((p,7)=9,P))/2+0LC 

S+(0.42<19)/Q 

A+ 2.50662823884 ~18.6150006252 41.39119773534 ~25.44106049637 
B+ ~8.4735109309 23.08336743743 ~21.06224101826 3.13082909833 
TeTx (C(T*2 )°.%0,13)+.%A F140 (7x2 owe LU dD +L XB 
Z0(0.422/Q)/1p,Q)+«7 

+(F=1)/0 


EXT:C+ ~2.78718931138 ~2.29796479134 4.85014127135 2.32121276858 


D+ 3.54388924762 1.63706781897 
S<(xS)«x((Re.*0,13)+.xC )#14+(((R«(]00.5-|S)*0.5)eo.% 1 2)4+.xD) 
Z0(0.42<1Q)/19QI)<S 


+0 

ERR:'ONE QR MORE P VALUES ARE QUT QF RANGE.' 
V 
V BETA (QO) 


DDDD BD 


V UtA BETA X;Y;W3N;OD;EV;Z31 


STATISTICS, 1978, P 325. BECAUSE OF THE RANGE OF 1, +/AS255. SEEMS IQ 


— oe ee ee ee ee eee oe os a —< === 


_———_ =e = —_— — oe oe —_— — — oe ot et ee ee ee oe 


Y+Xs(AC1)]#+/A) 

U+(p,X)p0 

N«744/(Cf /A)>(2*x14),10*110 

+((+/Y)=0)/FLIP 

WeY/X+«,X 

OD«Wo.x((1N)xAL2]-1N)#x/(N,2)pAC1]+12xI<Nn 

EV+-We.x(x/((2,N)p (AL1)4+0,1N-1), (+/A)4+0,1N-1))#x/(N,2)pAL1)]+0,1¢2*xN-Z+1 ) 


DPsZeitev ite +Ove si eZ 


+((I+«I-1)>0)/L 
ULY/1pUJ«(#Z)« (AL1) 7144/0 )x (W*Al1] )x(1-W)*A2) 
+((+/Y )=9X)/0 


FLIP:A«oA 


pe al Gad 86 
OD+We x (CiN)xAC2)~-1N)tx/(N, 2 )pAlL1) +12xI4Nn 
EV+-Wo.x(x#((2,N)p(AC1)40,1N-1),(+/A)4+0,1N-1) )#x/(N, 2 )pAl1)+0,1(2*xN-Z1 ) 


C21) £21:2Z<«1+EVC372141+0D0(0;72) +2 


(22) 
C23] 


+((T«I-1)>0)/L1 
UC CY )/19UJ)41-(4#Z)x(AL1) $714+4+/A )x (WxAL1) )x(1-W)*«AL2) 
v 
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